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ABSTRACT 



Q\ 

The role played by velocity fields in the fragmentation of a cold medium and in the 
£h formation of protostars is studied. 

I ^ The velocity field is modeled with a compressible turbulent flow. A supersonic turbulent 

velocity field can fragment the medium into clumps of mass smaller than a local Jeans' 
mass, and therefore stabilize the medium against the formation of protostars. Based 
on this idea, the protostar formation efficiency and the protostar mass distribution are 

£N| determined as functions of the following ambient parameters: average density no, aver- 

age temperature To, r.m.s. turbulent velocity a v o (or its Mach number Ait), postshock 
cooling time (e.g. chemistry). 
The main results are: 

o 

(i) the protostars mass distribution and its dependence on the ambient parameters are 
ON quantified; 

(ii) the characteristic protostar mass is Mj c \ oc n 1 ^ 2 Tq(t~ 1 q ; 

Oh (iii) the protostar formation efficiency e is higher for larger mean density, larger mean 

q temperature, lower velocity dispersion on a given scale and longer postshock cooling 

time (e.g. lower metallicity): e oc 1 ^Tq~ 1 (t v ^ ^ where f3 > 1 is the 

X/l exponent of the clump mass distribution; 

(iv) the efficiency is quite sensitive to the ambient parameters and therefore to the dy- 
namical evolution of the star forming system. 

Key words: stars: formation - turbulence 



1 INTRODUCTION 

1.1 Star formation in galaxies 

The observational properties of galaxies, upon which cosmo- 
logical models rely, cannot be interpreted without a knowl- 
edge of their star formation history. Examples of such prop- 
erties are integrated colors, surface brightness and chemical 
abundances. 

Also on purely theoretical grounds star formation is an 
essential part of the theory of galaxies formation, since the 
dissipative collapse of a protogalaxy is eventually arrested 
by its fragmentation into stars. Star formation is therefore 
a fundamental cosmological problem. 

The study of the evolution of a single protostar, while 



being a necessary part of the theory of star formation, is not 
the key point in the context of galaxy formation. Rather it 
is the statistical description of the fragmentation of a large 
scale environment with given initial and boundary condi- 
tions. Because of the complex network of interacting physi- 
cal phenomena in (proto)galaxies large scale protostars for- 
mation is certainly a chaotic process, so that its description 
must be statistical. 

Protostar formation in galaxies is normally described 
with simple parametrizations of the star formation efficiency 
such as Schmidt's law (Schmidt 1959; see Kennicutt 1989) 
or more complex (Wyse and Silk 1989; Dopita and Ryder 
1994), and often based on the idea of self regulated star 
formation (McKee and Lin 1988) and disk instabilities (see 
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the review by Elmegreen 1993). Feedback arguments are also 
often used (e.g. Whitworth 1979, Lin and Murray 1992). 

The problem is that of predicting protostar formation 
efficiency (or rates) and protostellar masses in a large scale 
environment with given dynamics and chemistry and which 
is initially devoid of stars. Without a statistical description 
of the formation of protostars our knowledge on protostar 
evolution will not have the desired impact upon the theory 
of galaxy formation and evolution. 

The present work is an attempt to predict protostar 
formation efficiency and protostellar masses as functions of 
ambient properties such as mean density, mean tempera- 
ture, turbulent velocity dispersion on a given scale, chem- 
istry (postshock cooling time) by focusing only on one aspect 
of the process of star formation, namely the hydrodynamics 
of the fragmentation of a cold medium when high Reynolds 
and Mach numbers flows are present. 



1.2 Supersonic turbulence and statistical 
description of protostars formation 

Many astrophysical environments possess highly supersonic 
flows with formation of shock waves. 

When the Reynolds number of the flow, that measures 
the ratio of the inertial term, v • Vv, relative to the viscous 
term, i/V 2 v, in the Navier-Stokes equation, is very large an 
important instability occurs: a regular flow (normally re- 
ferred to as laminar) is not possible, because the slightest 
disturbance will initiate the transfer of energy between dif- 
ferent scales in the flow, and the molecular viscosity cannot 
compete with this energy transfer due to the inertial term 
(see Batchelor 1967). The energy transfer occurs in the form 
of energy cascades, normally from the large scales to the 
small ones. 

The manifestation of this instability is generically called 
turbulence (see Batchelor 1967, Landau and Lifshitz 1987, 
Frisch and Orszag 1990). 

The importance of the onset of turbulence is that it al- 
lows for the use of equilibrium laws in the statistical descrip- 
tion of the flow on small scales (see Batchelor 1953, Monin 
and Yaglom 1971, Tennekes and Lumley 1972). Therefore 
while some astrophysical flows are certainly chaotic in the 
sense of not being predictable, others are chaotic and yet 
predictable in a statistical sense because they are unstable 
to turbulence. 

If this is the case for some astrophysical supersonic flows 
it is possible that the formation of shocks in them can be 
described as the effect of turbulence and therefore predicted 
in a statistical sense. Note that it was recognized long ago 
that the Reynolds number in the interstellar matter on the 
scale of lOOpc has a very high value of the order of 10 6 (von 
Weizsker 1951). 

The problem is particularly interesting in the star form- 
ing environments. Understanding protostar formation is 
mainly a matter of predicting the local occurrence of grav- 
itational instabilities in large regions of the interstellar gas, 
not necessarily unstable as a whole. The simplest way to 
do this is to assume that the gas fragments into protostars 
whose mass equals the Jeans' mass in the large region which 
is supposed to be homogeneous and with little velocity fluc- 
tuations (Jeans 1902, Hoyle 1953, Fowler and Hoyle 1963). 



This approach is not even a first approximation if the 
interstellar gas is organized in supersonic chaotic flows be- 
cause of the large velocity and density gradients. 

Nevertheless if the motions are in the regime of fully de- 
veloped turbulence it could be possible, as mentioned above, 
to predict the formation of shocks using the statistical the- 
ory of turbulence and apply simple instability arguments to 
the shocked regions. 

In this way the protostar formation will be shown to 
depend not only on the mean density and temperature of 
the environment, but also on the turbulent velocity disper- 
sion responsible for the compressions in the gas, and on the 
chemistry responsible for the postshock cooling. Clearly the 
more that is known about the turbulent model flow, the 
more can be inferred about the shock distribution and the 
protostar formation. 

This simple way of describing protostar formation is 
certainly not comprehensive of all the physics involved in 
the problem of large scale star formation (especially because 
magnetic fields will not be considered in the following). Nev- 
ertheless it is interesting to study the role of hydrodynamical 
supersonic turbulence in the formation of protostars, at least 
in environments where the magnetic field is known not to 
play a dominant part. 

Larson (1981) suggested that present day star forma- 
tion in the Galaxy could occur in shocks generated by the 
supersonic turbulence that could be the reason for the width 
of the molecular lines emitted by the interstellar gas. 

Hunter (1979) and Hunter and Fleck (1982) have stud- 
ied the influence of the velocity field upon star formation, 
but without taking into account the formation of shocks. 

Leorat et al. (1990) have shown, with two-dimensional 
numerical simulations, that compressible turbulence can af- 
fect the gravitational instability in the ISM. 

Elmegreen (1993) has studied the gravitational instabil- 
ity of shocked regions of the ISM and shown how the clumpy 
structure of molecular clouds and their star formation mech- 
anism can be controlled by the turbulent Mach number. 



1.3 Overview of the paper 

Section 2 contains a discussion of the physical assumptions 
used throughout the work. 

The fragmentation controlled by the turbulent velocity 
field is described in terms of postshock high density clumps 



* Chandrasekhar (1958) suggested a correction of the Jeans' mass 
due to turbulent pressure; Arny (1971) used the idea of turbulent 
pressure to predict the protostar mass function; Bonazzola et al. 
(1987) studied the influence of the turbulent spectrum on the 
gravitational instability, restudying the effect of turbulent pres- 
sure. That turbulence acts as a pressure source against gravity has 
been shown by Bonazzola et al. (1992) (see also the simulations 
by Leorat et al. 1990). 

Nevertheless, in order to be the dominant source of pressure 
the turbulent motions must be supersonic (only in this case their 
role is significant in the cited works); if this is the case and if 
the medium has a short cooling time (in the sense specified in 
the paper) strong shock formation leading to very high density 
contrasts is likely to be the main fragmentation mechanism due 
to turbulence, and turbulent pressure is of secondary importance. 
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in section 3, where the properties and stability of the clumps 
are studied. 

The concept of turbulent fragmentation is applied to 
protostar formation in section 4; in particular formulae for 
the mass distribution and the formation efficiency of proto- 
stars are given. 

Section 5 is a discussion of the role of turbulent frag- 
mentation in controlling the formation of protostars. 

The work is summarized in section 6. 

A list of the symbols used throughout the paper is given 
in the appendix. 



2 PHYSICAL ASSUMPTIONS 

A number of physical assumptions will be used in order to 
study the fragmentation controlled by a turbulent velocity 
field. 

The first is the exclusion of the magnetic field. The 
observational data on the intensity of the magnetic field 
in dense clouds (Heiles and Stevens f986; Crutcher, Kazes 
and Troland f987; Kazes et al. f988; Crutcher et al. f993) 
would not justify this assumption if the present Galactic 
ISM were the subject of the study; nevertheless it is im- 
portant to understand the role played by purely hydrody- 
namical processes in the fragmentation of a cold medium. 
The whole hydromagnetic treatment of the problem is under 
study through numerical simulations and will be presented 
in a later paper (Nordlund and Padoan, in preparation). 

It will be shown that if a protostar is identified with one 
Jeans' mass in the fragmented gas, the protostars mass dis- 
tribution and formation rates can be determined once the 
energy spectrum of turbulence and the mass and density 
distribution of postshock clumps are specified. Note that 
the essential ingredient is the energy spectrum of turbu- 
lence, since this fixes the typical density and mass of the 
clumps and the typical protostellar mass. It is therefore as- 
sumed that the ISM flows are in the regime of fully developed 
turbulence, as their huge Reynolds number indicates (Von 
Weiszker 1951; Scalo 1987), so that a power law velocity 
spectrum can be assumed. 



2.1 The picture 

A complex system of isothermal shock waves arises in the 
ISM due to the presence of a turbulent and highly supersonic 
velocity field. After one dynamical time of the largest scale 
the ISM matter distribution is very clumpy, with density 
contrasts of the order of the mean squared Mach number of 
the velocity fluctuations. For the sake of definitiveness it will 
be assumed in the following that the mass distribution can 
be described as an ensemble of clumps that are formed after 
shock waves. A typical mass scale of turbulent fragmentation 
will be called "typical clump" and a clump mass distribution 
will be defined. 

Some clumps happen to be dense and massive enough 
to undergo gravitational instability and fragment into pro- 
tostars. Clearly this way of forming protostars is highly de- 
pendent on the density and on the mass of the clumps. More 
precisely, since in typical ISM conditions most clumps will be 
stable against gravitational collapse, it is the intermittency 
in the statistics of the density field (shaped by supersonic 



turbulence) that determines what fraction of the total mass 
becomes unstable and is converted into protostars. 

2.2 Turbulent spectrum 

Supersonic velocity dispersions a v are observed on scale r, 
according to the Larson's relation (Larson 1981): 



lkms 



lpc 



(1) 



valid in the range of scales O.f — I OOpc (Larson 1981). By in- 
cluding HI regions this relation is valid up to fOOOpc (Larson 
1979). The correlation between line width and linear size in 
the molecular ISM has been confirmed (with slightly differ- 
ent exponents and normalization factors) by a number of au- 
thors (Leung, Kutner and Mead 1982; Myers 1983; Quiroga 
1983; Sanders, Scoville and Solomon 1985; Arquilla and 
Goldsmith 1985; Dame et al. 1986; Falgarone and Perault 
1987; Fuller and Myers 1992). 

Despite the numerous interpretations of this relation 
in terms of a Kolmogorov turbulent cascade (Kolmogorov 
1941), it should be remembered that Larson's relation is 
not a direct measure of the velocity power spectrum, so that 
the observed exponent, 0.4 — 0.5, could arise from different 
mechanisms rather than from the presence of a turbulent 
spectrum such as the one expected from models of turbu- 
lence. 

On the other hand a number of authors have tried to 
measure autocorrelation and structure functions of ISM ve- 
locity fields directly (Scalo 1984; Dickman and Kleiner 1985; 
Kleiner and Dickman 1985; Kleiner and Dickman 1987; Hob- 
son 1992, Miesch and Bally 1993). No definitive conclusion 
as to which model of turbulence is most appropriate to de- 
scribe the ISM has been reached. 

When modeling the ISM turbulence one has therefore 
to rely on the statistical description of numerical turbulence. 
Direct numerical simulations of highly supersonic turbulence 
are still very difficult, but there are good indications that an 
energy spectrum: 



E(k) oc k~ 



(2) 



is to be expected for the compressional modes (Passot, Pou- 
quet and Woodward 1988; Porter, Pouquet and Woodward 
1992) that are responsible for the fragmentation (formation 
of shocks). Such a spectrum is also to be expected for a flow 
made up of discontinuities (Passot and Pouquet 1987), and 
it is known to arise from Burgers' equation (Burgers 1940, 
1974; Gotoh and Kraichnan 1993). 

Given these indications, the power spectrum (2) will be 
assumed. 

In equation (2) E(k) is the integration, over a spherical 
shell of radius k (modulus of the wave vector), of the trace 
$ 88 (k) of the spectrum tensor $ 8J (k) which is the Fourier 
transform of the correlation tensor i? 8J (r): 



Rzj(r) = m,(x, t)Mj(x + r, t) 



(3) 



where u t is the flow velocity. 

The trace $ 88 (k) represents the kinetic energy per unit 
mass at wave vector k, as can be seen from: 



Ra(o) 



u,u, 



$„(k)dk 



(4) 
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Thus one can write: 

f CO 

n l n l oc / E(k')dk' oc k -1 oc r 

Jk 



(5) 



Therefore if the spectrum (2) is used, the r.m.s. velocity 
cr v (r) on the scale r (or the velocity of "eddies" of size r) is: 



/ \ i/ 2 
cr v (r) oc r 



(6) 



Notice that this relation means that the dynamical 
time, r/a v (r), increases with the linear scale according to: 

Tdyn(r) OC f 1/2 (7) 



2.3 Clumps density and mass distributions 

The shock waves in the flow result in very non-Gaussian 
statistics for the velocity gradient. The statistics of diver- 
gence of the velocity field has been shown to be extremely 
intermittent in numerical simulations of compressible tur- 
bulence (Lee, Lele and Moin 1991). The density distribu- 
tion is thus expected to be non-Gaussian; Vazquez-Semadeni 
(1994) has found that the density field in two-dimensional 
simulations of compressible turbulence is lognormally dis- 
tributed; this is confirmed in preliminary results of three- 
dimensional numerical simulations of highly supersonic tur- 
bulence (Nordlund and Padoan, in preparation). 

Therefore the density distribution of the clumps into 
which the postshock gas is fragmented is assumed to be log- 
normal: 



P(x)dx 



(27T) 1 / 2 



exp 



- — ln 2 x \ dx 
2 



(8) 



where x = n/n r j, with n number density and n r j the density 
of the typical clumps (defined by equation (10) that depends 
on the normalization of the turbulent spectrum (see next 
section). 

It will also be assumed that the density and mass distri- 
butions are statistically independent so that the joint prob- 
ability can be written as a product of the two distributions 
(this will simplify the derivation of the protostars mass dis- 
tribution). This is justified by the fact that the density is 
mainly determined by intense shocks on large scales while 
the fragmentation into small clumps is controlled by the less 
intense shocks on small scales. 

The clumps into which the dense postshock gas is frag- 
mented are characterized by a typical mass, M r j (see equa- 
tion (15), that depends on the turbulent spectrum (see next 
section); nevertheless it is clear that larger clumps can exist, 
even if less numerous. For simplicity the mass distribution 
is assumed to be a power law: 



Pci(m c i)dm c i = (1 — fi)mj' dm c i 



(9) 



for rrid > 1, where m c i is the mass in units of the typical 
clump mass, M r j, defined in the next section, and fi > 1. 

Note that this mass distribution (as all the following 
ones) is defined in such a way that the number of objects 
per mass interval is N(m)dm oc P ^™' 1 dm, that is P(m) is 
proportional to the number of objects per logarithmic mass 
interval. 

Several observations yield the value fi PS 0.6, for both 
giant molecular clouds and clouds cores. The subsection 5.1 



explains why there is no contradiction between the observa- 
tions and the assumed clump mass distribution. 



3 TURBULENT FRAGMENTATION 
3.1 Clumps properties 

Given the turbulent velocity (6), it is possible to define the 
transition scale rtr below which the velocity becomes sub- 
sonic: 



M(r tr ) = 1 



(10) 



where M(r) is the Mach number of the turbulent velocity 
dispersion on the scale r, relative to the isothermal sound 
speed of the medium. 

The scale rtr is the typical scale of the fragmentation 
into clumps due to a turbulent velocity field. In a cloud of 
mass Mdoud it corresponds to a mass given by: 

Mdoud 



Mtr 



M* 



(11) 



where Mt is the Mach number of the cloud velocity disper- 
sion; Mtr is between 0.01 and IOMq in typical giant molec- 
ular clouds of the galactic disc. This is an interesting range 
of masses for protostar formation. 

If the postshock cooling time is supposed to be shorter 
than the dynamical time of the shocks (see below), the post- 
shock density on the scale rtr is estimated by the isothermal 
jump condition: 



— ~ Mt 

Pi 



cs 



(12) 



where p\ and p2 are the gas volume density before and after 
the shock, c s is the isothermal sound speed and Mt is the 
Mach number of the ambient r.m.s. turbulent velocity, <J v ,o, 
on the scale Lq. 

The reason why this large scale Mach number is used 
is that the stability of the shocked structures formed after 
one dynamical time of the large scales will be later studied 
(this will be used to define the efficiency of the formation of 
protostars). In fact most of the density contrast relative to 
the mean density is produced by large scale shocks after one 
dynamical time of the large scales. Shocks on small scales 
are less intense (see equation (6)) and their role is mainly 
that of fragmenting the gas into smaller subunits, since their 
contribution to the density enhancement over the mean am- 
bient density is irrelevant. 

The relations (6), (10) and (12) determine the mass and 
the density of the typical postshock clumps into which the 
medium is fragmented. The typical mass M r j is: 



,t 4-7T 3 

Mcl = —Por tr 



(13) 



where po is the ambient average density and the transition 
scale defined by (10) is: 



r t 



L 



(14) 



lOOpc 

where Lo is the linear size of the region whose r.m.s. turbu- 
lent velocity is a Vi o (e.g. the size of the energy containing 
turbulent eddies). 
Therefore: 
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Figure 1. The typical mass of the postshock clumps, M c i , and the Jeans' mass in the clumps, Mj >c i , are plotted versus the ambient r.m.s. 
turbulent velocity (a), the mean density (b) and the mean temperature (c). In (a) no = 100cm -3 and To = 10°A , in (b) <r V) o = 2km/ s 
and To = 10°A , in (c) <r V) o = 2km/ s and no = 100cm -3 have been used. In all three plots To = lOOpc. The mass is expressed in solar 
masses. 



M « « 1O ^0(T^r^)(^77)" 6 (TSf) 3 (T^) 3 

«i°^0(i^^)(^r 6 (i^) 3 ( 15 ) 

where no is the number density corresponding to po. 

From the relation (12) the number density in the clump 

is: 

n r j « n M 2 t (16) 

If the postshock cooling time, r coo i, is longer than the 
shock dynamical time, Td yn , (the time required by the shock 
to traverse the clump's mass) no postshock structure of high 
density contrast of mass comparable to the one traversed 
by the shock can be formed. The radiative postshock zone 
would extend through the whole shocked gas, and the gas 
would simply disperse, since the energy of compression is 
expected to be larger than the gravitational energy on this 
scales (see Larson 1988). Because of the turbulent energy 
spectrum k~ 2 (equation (2)) this timing condition for the 
shocks to be radiative can easily be translated into a mini- 
mum mass, M ra d, for postshock clump formation: 

M rad « Mo (IzzL) (17) 

\Tdyn,0 J 

where Mo is the mass contained in the large scale with dy- 
namical time Tdynfi- 

Since for a large system r coo i < Tdynfi (see equation 
(7)) and because of the large exponent in the relation (17) 
M r ad <C Mo is generally satisfied. 

The typical mass of the fragments produced by the tur- 
bulent velocity field is therefore max(M c i, M ra d)- For a so- 
lar metallicity the cooling time in molecular clouds is short 
enough that M ra d < M r j. The low metallicity regime with 
M r ad > M r j is of great importance to describe population 



II star formation. There are a number of interesting issues 
relating to this regime and this will be dealt with in a later 
paper. 

3.2 Clumps stability 

3.2.1 Turbulent pressure 

Molecular clouds have too low temperatures (10 — 50K) 
to be supported by their thermal pressure (Larson 1981; 
Bonazzola et al. 1987, 1992; Leorat, Passot and Pouquet 
1990). The ISM turbulence has been invoked as an al- 
ternative to the magnetic field (see Mouschovias 1987) as 
the pressure source against the gravitational collapse of 
molecular clouds and different mechanisms of energy injec- 
tion into the turbulent cascade have been proposed (von 
Weizsaecker 1951; Fleck 1980,1981; Norman and Silk 1980; 
Scalo and Pumphrey 1982; Franco and Cox 1983; Henrik- 
sen and Turner 1984; McCray and Kafatos 1987). A pos- 
sible solution to the problem of the short dissipation time 
of the turbulent motions compared to the estimated life- 
time of giant molecular clouds (Blitz and Shu 1980), is that 
the dissipation of turbulence is reduced as a consequence 
of the fragmentation of the clouds into clumps (Scalo and 
Pumphrey 1982). 

Turbulent pressure could be important also inside the 
postshock clumps into which the medium is fragmented by 
supersonic turbulence (see subsection 3.1), but it is probably 
not the case because of the following arguments. 

The density contrast between the clumps and the sur- 
rounding medium and the resultant low volume filling frac- 
tion of the clumps is likely to interrupt the hydrodynamic 
turbulent cascade of energy from the large scales. Moreover 
the fragmentation inside the clumps cannot be invoked to 
make the dissipation time longer, since for the definition 
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Figure 2. Protostars mass distribution from equation (18). The function is plotted for different values of the ambient density (in cm -3 ) 
(a), turbulent velocity dispersion (in km/ s) (b) and temperature (in K) (c). In (a) To = 10A' and <J V) o = 4km/ s; in (b) no = 200cm -3 
and To = 10A ; in (c) no = 200cm -3 and <J V) o = 4km/ s. The mass is in solar masses. The clumps mass distribution has been modeled 
with a power law (see text) with exponent f3 = 1.5 



(10) of the clumps no further fragmentation will occur in- 
side them. Thus the turbulence inside the clumps is probably 
dissipated in a short time (of the order of, or less than, the 
dynamical time in the clumps). Finally turbulent pressure 
is not dominant over thermal pressure (as is obvious from 
the definition (10)), despite some increase of turbulent en- 
ergy in the postshock gas is expected (see Rotman 1991, 
Jacquin, Cambon and Blin 1993). For these reasons the in- 
ternal turbulent pressure in the clumps will not be used in 
the definition of the Jeans' mass inside the clumps. 



3.2.2 Jeans' mass in the clumps 
The Jeans' mass inside the clumps is: 



becomes larger than the largest stable mass, Mj tC i, and 
therefore the clumps are unstable, for decreasing ambient 
r.m.s. turbulent velocity and for increasing mean density and 
temperature. The increase in a Vi o produces denser clumps 
(higher Mt), but since their typical mass is also reduced, 
they are more stable. 

If the postshock cooling time is such that M ra d > M r j 
(see eq. (17)), the mass to be compared with Mj tC i is M ra d 
and not M r j. In this regime a long postshock cooling time 
means more massive and more unstable clumps than in the 
case of a shorter cooling time. Inefficient cooling acts against 
the fragmentation by the velocity field that otherwise tends 
to stabilize the medium against the gravitational instability 
by producing clumps smaller than the local Jeans' mass. 



Mj, cl ^UMrAjhf 1 ^^^)- 112 ^)- 112 

« 12M ( T X f ) 2 ( TIJI ^ FT ) -1 / 2 (^ 7 ) -1 (18) 

where equation (16) has been used in the second expression 
and equation (12) (the definition of the Mach number) in 
the third, with the result of changing the temperature de- 
pendence and of introducing the ambient velocity dispersion. 
This Jeans' mass decreases with increasing ambient velocity 
dispersion (denser clumps are produced according to equa- 
tion (16)) contrary to the case of internal turbulent pressure 
(see references in the footnote in section 1). This is in fact 
equivalent to the Jeans' mass at constant external pressure, 
if turbulent ram pressure is considered. 

To study the stability of the typical clumps formed by 
the supersonic velocity field, equation (18) has to be com- 
pared with equation (15). This is illustrated in figs. 1. 

As shown in figs.l the mass of the typical clumps M r j 



4 PROTOSTAR FORMATION 

4.1 Protostar mass distribution 

It is assumed that any unstable clump fragments into proto- 
stars whose mass equals the Jeans' mass inside the clump. ^ 
Since a density and a mass distributions of clumps are as- 
sumed to exist, with typical values n r j and M r j determined 
by the turbulent spectrum, a protostar mass distribution 
will also exist around the value Mj tC i- 

If n r j in equation (18) is changed into n, that is the local 
number density and obeys the statistics (8), equation (18) 

t Thin slabs are more stable than the simple Jeans' mass pre- 
dicts (see for example Gilden 1984); nevertheless in this work it 
is assumed that the velocity field fragments the medium down to 
the small scales where the typical shocks are not very intense and 
therefore the postshock gas cannot be described as a thin slab. 
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Figure 3. The protostar formation efficiency per ambient dynamical time is plotted versus the ambient density expressed in cm -3 
(a), the turbulent velocity dispersion in km/s (b) and the ambient temperature in K (c). In (a) <J V; o = 4km/ s and To = 10A ; in (b) 
no = 100cm -3 and To = 10A ; in (c) no = 100cm -3 and <J V) o = 4km/ s. To = lOOpc and (3 = 2.0 in all three plots. The curves are power 
laws (see text). 



gives the mass Mj of the protostars formed in the clumps 
of density n. 

The protostars mass distribution P(m), with m = 
Mj/Mj t ci, is therefore determined by: 

P( m ) = P{x)^- (19) 

(x = n/n r j) where P(x) is from equation (8) and dx j 'dm 
comes from equation (18) after the substitutions Mj tC i = Mj 
and n r j = n. The protostar mass distribution is thus: 

2m — 1 

P(m)dm = - — — —r- exp ( — 2 In m)dm (20) 

(27T) 1 ' 2 

If the mass distribution of clumps is P c i(m c i), the pro- 
tostar mass distribution becomes: 

/CO 
P c i(tn c i)dm c i (21) 

which means that protostars of mass Mj can be formed only 
in clumps of mass > Mj = mMj iC i/M r j. 

The factorization into the product of the protostar mass 
distribution and the integral of the clumps mass distribution 
is allowed by the assumption of statistical independence of 
the density and mass distributions of the clumps (see sub- 
section 2.2). 

Using the clumps mass distribution (9), the protostar 
mass distribution, equation (21), is further transformed into: 

P(m)dm = - — — — - exp( — 2 In m)dm, 

(27T) 1 ' 2 

ifm < M c i/Mj tC i; 

u( \a 2m ~' 3 i Mcl ^-i ( 11 2 \A 

P{m)dm = —^—{jj—y exp(-2 In m)dm, 



ifm> M c i/Mj,d (22) 

Several examples of protostar mass distributions are 
plotted in fig.s 2, for different values of the ambient param- 
eters. 

In each of fig.s 2 three protostar mass distributions are 
plotted for different values of no, <Jv,o and To respectively. 
The function plotted is equation (22) with fi = 1.5. 

The protostar mass distribution has a single maximum 
corresponding approximately to 0.5Afj |C i, and goes to zero 
for both very large and very small protostellar masses. 

It has to be emphasised that the characteristic mass of a 
protostar given by equation (18) differs from the traditional 
Jeans' mass in the temperature dependence and because it 
depends also on the ambient velocity dispersion cr Vi o- 

Mj iC i oc n^ l/ ' ''To <t~o (23) 

This characteristic protostellar mass decreases with in- 
creasing velocity dispersion (because denser clumps are pro- 
duced according to (16)), and is equivalent to the Jeans' 
mass at constant external pressure, if turbulent ram pres- 
sure is considered. 

If the postshock cooling time is such that M ra d > M r j 
(see equation (17)) the protostar mass distribution is the 
same as in equation (22), but with the substitution M r j = 
M ra d, because no high density contrast bound clump of mass 
< M r ad can be formed, as explained in the end of the sub- 
section 3.1. 

4.2 Protostar formation efficiency 

The protostar formation efficiency is simply determined by 
the mass fraction of unstable clumps relative to the total 
mass. The discussion about the postshock clumps stability 
presented in section 3 is therefore equivalent to the present 
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discussion of the efficiency of protostar formation. The sce- 
nario is now complicated by the introduction of a mass dis- 
tribution for the clumps and of a distribution for the Jeans' 
mass inside them (that are supposed to be statistically in- 
dependent). 

The mass distribution (21) expresses the fraction of the 
total mass contained in clumps whose Jeans' mass relative 
to Mj t ci is between m and m + dm and whose mass is > m. 
Therefore the integration of (21) over m gives the fraction 
of total mass contained in unstable clumps (or equivalently 
in protostars): 

edyn = * = / P{m)dm (24) 

Mtot J 

The compression ratio cx M 2 t , used to determine the 
clumps density (see section 3), is experienced by the whole 
environment after one dynamical time of the large scale Lo : 

Tdynfi = — (25) 

<?vfi 

where a Vi o is the r.m.s. velocity on the scale Lo. Therefore a 
time equal to Tdynfi is necessary to convert into protostars 
as much mass as expressed by equation (24). 

For this reason the efficiency, defined as the protostar 
formation rate per unit mass, is: 



In figs. 3 the dependence of the efficiency on ambient pa- 
rameters is shown: the efficiency per dynamical time edyn, 
equation (24), is plotted versus the mean density, the am- 
bient r.m.s. turbulent velocity and the mean temperature 
respectively. The clumps mass distribution is modeled with 
a power law with exponent fi = 1.5 and M ra d < M c i is 
supposed to be valid. 

The efficiency grows with no and To and with decreasing 
<T V fi, as the ratio M r j/Mj iC i (see subsection 3.2). 

The curves in fig. 3 follow the power laws: 

from edyn ~ 0.01 up to values as high as edyn ~ 0.9. 

This can be also seen qualitatively since the factor that 
multiplies the protostar mass distribution equation (22) for 
m > M c i/Mj tC i, (Md/Mj^i) 13 ' 1 , has exactly the same de- 
pendence on the ambient parameters as edyn expressed in 
equation (27). 

It is interesting that while it is possible that the effi- 
ciency per dynamical time PS 3% at To = 10i<f , the efficiency 
in the same environment warmed up to To = 30 K increases 
to PS 9% (for <T V fi = Akmjs on the scale To = lOOpc and 
no = 100cm -3 ). This fact seems to indicate that wherever 
the mode of protostar formation here discussed is relevant, 
the protostar formation rate can be increased by the warm- 
ing up of the ambient medium due to the first generation 
of stars (positive feedback). The increased temperature will 
also increase the protostars mass (fig. 2.c) with the effect 
of an even more efficient warming of the environment; this 
makes the feedback stronger (and the characteristic proto- 
star mass increase with time). 

A negative feedback could instead result from the injec- 
tion of mechanical energy due to the products of protostar 
formation, if this could increase the a Vi o (see fig. 3.b). 



If the ambient chemistry is such that M ra d > M c i (see 
equation (17)) M ra d should be used instead of M r j for esti- 
mating the efficiency; if M ra d/M c i 3> 1 the increase in the 
efficiency could be drastic. 



5 DISCUSSION 

It is a common assumption in theories of galaxy formation 
that, when the cooling time of the gas is shorter than the free 
fall time of the protogalaxy, fragmentation into stars occurs 
(Hoyle 1953; White and Rees 1978). This notion is based on 
linear gravitational instability and its applicability to a flow 
with certainly nonlinear velocity and density fluctuations is 
dubious. 

Since the flow is turbulent and supersonic the study 
of the fragmentation of the protogalactic gas when cool 
and with short cooling time is a problem of hydrodynamics 
rather than of gravitation, as far as the baryonic component 
of the galaxy and small scales are concerned. 

The idea discussed in this work is that the fragmenta- 
tion of a cold medium is first controlled by the supersonic 
turbulent velocity field, while gravity takes over only inter- 
mittently. The large scale velocity field is determined by the 
gravitational potential, but on scales much smaller than the 
size of the galaxy the flow arises from hydrodynamical in- 
teractions inside the large scale flow; these interactions ap- 
pear as energy cascades (the Reynolds number is very large) 
that are able to feed turbulent velocity fluctuations on small 
scales. 

The velocity fluctuations are supersonic so that the 
fragmentation of the density field has little to do with grav- 
itational forces, at this stage. In this way a very clumpy 
density distribution arises out of a complex system of shock 
waves; it is on such a distribution that gravity applies, forc- 
ing the largest or densest clumps to fragment into protostars. 

If the first stage of the fragmentation based on turbu- 
lent supersonic cascades produces only very small clumps, 
as in the cold Galactic ISM (fig. 4), most of these will be sta- 
ble against the gravitational collapse. In this case it is the 
intermittency in the statistics of the density field (shaped 
by supersonic turbulence) that determines what fraction of 
the total mass becomes unstable and is converted into pro- 
tostars. 

Therefore turbulence could be the reason why present 
day star formation is observed to be a slow process, capable 
of converting only a few percent of the gaseous mass into 
protostars (Duerr, Imhoff and Lada 1982; Myers et al. 1986; 
Mooney and Solomon 1988) in a period of about 10 7 years 
(Blitz and Shu 1980) after which the stellar feedback de- 
stroys the star formation site (Herbig 1962; Field and Saslaw 
1965; Elmegreen and Lada 1977; Whitworth 1979; Lada and 
Gautier 1982). In fact the formula (27) for the protostar for- 
mation efficiency derived here predicts the conversion into 
protostars of only a few percent of the total mass after 10 r 
years, if applied to Galactic molecular clouds. Moreover the 
predicted mass distribution can be close to the observed 
Miller-Scalo mass distribution (Miller and Scalo 1979). 

The main achievement of the present work is that of 
showing that supersonic turbulence makes gravitational frag- 
mentation and thus star formation inefficient with the result 
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Figure 4. The turbulent fragmentation mass scale is compared 
with the Jeans' mass in the fragmented gas in the plane density- 
velocity dispersion. M c i = Mj c i along the line; the region above 
the line produces turbulent fragments that are smaller than their 
Jeans' mass. The position of giant molecular clouds and their 
cores is shown to be in the stable part of the plane, meaning 
that gravitational fragmentation of the cold ISM is hindered by 
supersonic turbulence. The cores position is at the same velocity 
as the parent GMC because the velocity is referred to a fixed large 
scale (this is a consequence of Larson's relation). 

that the protostar formation rate increases with increasing 
average temperature and decreasing velocity dispersion. 

Since the role of the magnetic field has not been consid- 
ered one should be cautious when trying to predict present 
day star formation using the formulae here derived; never- 
theless the present work shows that the turbulent fragmen- 
tation is certainly a powerful mechanism to slow down the 
process of star formation and even to give the stellar mass 
distribution its initial shape. 

5.1 ISM dumpiness 

The present model relies on the basic idea that the ISM, 
due to ubiquitous supersonic turbulent flows, is fragmented 
down to the scale defined by equation (10), that marks the 
transition to subsonic motions. Therefore most of the mass is 
supposed to be contained in fragments of that typical scale, 
so that if a mass distribution is assumed its exponent has to 
be /3 > 1. 

Most analysis of molecular cloud maps have revealed 
the presence of a very clumpy structure, with clumps mass 
distributed on power laws with exponent (defined in subsec- 
tion 5.3) 0.0 < $ < 0.7 (Myers, Linke and Benson 1983; Blitz 
1987; Carr 1987; Loren 1989a; Stutzkie and Glisten 1990; 
Lada, Bally and Stark 1991; Nozawa et al. 1991; Langer, 
Wilson and Anderson 1993; Williams and Blitz 1993). 

This means that most of the detected CO mass is con- 
tained in the few largest clumps identified, with mass of 
the order of 100 — IOOOMq and internal velocity dispersion 
larger than lkm/ sec. These distributions are determined for 
clumps in the mass range 10 — IOOOMq, because of incom- 
plete sampling for smaller clumps, and always detect cores 
with supersonic velocity dispersion. 

Therefore almost none of the identified cores is a clump 
as defined in the present work (equation (10)), that is a frag- 
ment with nonsupersonic internal velocity dispersion: most 
of the observed cores must be internally fragmented. The 



predicted mass scale of turbulent fragmentation in the ISM 
is of the order of O.IMq, that is in fact two orders of magni- 
tude smaller than the lower limit of most mass distributions 
which have been reported in the litterature. 

Moreover it has been recognized that the ISM exhibits 
a hierarchical structure (Scalo 1985; Falgarone and Perault 
1987; Vazquez-Semadeni 1994). This fact alone is sufficient 
to prove that the mass distributions obtained from cloud 
maps are not intrinsic distributions of fundamental frag- 
ments, but rather the result of the selection of some step 
in the hierarchy as the consequence of a combination of ob- 
servational resolution, blending, projection and cores defini- 
tion. 

For these reasons the mass distributions obtained from 
the analysis of molecular cloud maps cannot be compared 
with the mass distribution used in the present work. 

On the other hand high resolution observations (Walm- 
sley et al. 1995; Velusamy, Kuiper and Langer, 1995; Langer 
et al., 1995) are proving that dark cloud cores are in fact 
fragmented down to very small scale into dense, unbound, 
young clumps with subsonic internal velocity dispersion, 
as expected from the present theoretical scenario. These 
clumps are transient structures not formed by gravity. 

High resolution observations of very high density trac- 
ers are not available yet for a statistical analysis of a large 
number of such small clumps whose mass distribution has 
never been observationally determined. 



5.2 Protostar formation 

A sub- Jeans fragment of shocked gas is in general a tran- 
sient structure in the ISM density field. Supersonic turbu- 
lence produces continuously such transient structures that 
are not able to collapse, as shown in fig. 4. This is a natural 
explanation, not invoking magnetic fields, of the fact that 
giant molecular clouds are not turning efficiently into stars 
in one free-fall time although they contain several Jeans' 
masses: supersonic turbulence makes gravitational instabil- 
ity behave as an intermittent process. 

The protostar mass distribution is mainly determined 
by the p.d.f. of the turbulent density field that is expected 
to be close to lognormal. This fact explains why the stellar 
mass function can be much steeper than the mass distribu- 
tion of cloud cores (see the previous subsection), although 
an alternative explanation is given in Zinnecker (1989). In 
any of the observed cores several stars can be formed with 
masses determined by the universal statistical properties of 
turbulence. 

Several complications can be included for a better de- 
scription of the formation of protostars and were not con- 
sidered in this work. 

First of all one Jeans' mass can fragment into a multiple 
stellar system instead of becoming a single protostar, and in 
this process only some fraction of the total mass will go into 
the final stars. 

Secondly, gravitationally bound groups of transient 
clumps can be formed. Clumps in such groups can coalesce 
and form an unstable core that will collapse. Therefore when 
the present model predict no star formation, clumps coales- 
cence could be the main star formation mechanism. Notice 
however that the collisional build up of an unstable core 
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must be the consequence of the dissipative contraction of a 
bound system of smaller clumps, and not of collisions in gen- 
eral, which are just one aspect of the supersonic turbulent 
flows that are here described. Since an appreciable fraction 
of identified cores seem to be unbound systems (for exam- 
ple in Blitz 1987; Carr 1987; Loren 1989b; Herbertz, Un- 
gerechts and Winnewisser 1991; Blitz 1993), the protostar 
formation via clump coalescence might be very inefficient. 
The largest cores are more likely to be gravitationally bound, 
but since their internal velocity dispersion is supersonic, tur- 
bulent fragmentation operates against protostar formation 
also inside these cores. 



no other mode will be important because the medium is 
expected to be continuously disrupted into very small frag- 
ments (compared to some local Jeans' mass). The only re- 
striction for the validity of the results here obtained is there- 
fore that the magnetic pressure has to be smaller than the 
thermal pressure. Since the cosmological seed magnetic field 
is expected to be < 10 -17 G (Lazarian 1992), it can be ar- 
gued that the present study is relevant to the problem of pri- 
mordial star formation that occurs at epochs when no large 
scale dynamo has yet been activated in galaxies (Moffat 
1978, Parker 1979, Krause et al. 1980, Ruzmaikin, Shukurov 
and Sokoloff 1988). 



6 SUMMARY AND CONCLUSIONS 

A supersonic turbulent flow has been used as a model for 
the velocity field in the cold ISM. 

A typical mass of postshock clumps is defined and the 
stability of these to the gravitational collapse is studied. 

The unstable postshock clumps are further fragmented 
into protostars by gravitational instability. 

A turbulent spectrum and a mass and density distribu- 
tion of the clumps are used to quantify the protostar mass 
spectrum and formation efficiency. It is shown how these 
depend on the ambient parameters. 

The main aspects of the present statistical description 
of protostar formation can be summarized in the following 
points: 

(i) the protostar mass distribution has a single maximum 
and depends on the ambient parameters. 

(ii) The characteristic protostar has a mass Mj tC i oc 
n 1 ^ 2 Tq Uy 1 , different from the traditional Jeans' mass 
and from the Jeans' mass with turbulent internal pres- 
sure, and equivalent to the Jeans' mass at constant ex- 
ternal turbulent pressure. 

(iii) The protostar formation efficiency is higher for larger 
mean density, larger temperature, lower turbulent ve- 
locity dispersion and, in the regime M ra d > M r j, 
longer postshock cooling time (e.g. very low metal- 
licity). The following power laws are followed: e cx 

4 (/3 " 1) To' 3 - 1 <o (/3 " 1) £o (/3 " 1 \ where /? > 1 is the ex- 
ponent of the clumps mass distribution. 

(iv) the efficiency is quite sensitive to the ambient param- 
eters and therefore to the dynamical evolution of the 
star forming system. 

The mode of protostar formation here discussed is rel- 
evant to environments with supersonic turbulent velocity 
fields, postshock cooling time shorter than the large scale 
dynamical time and magnetic pressure not larger than ther- 
mal pressure. Whenever these conditions are met it is likely 
that the study of highly supersonic turbulence will give con- 
siderable insight into the problem of large scale protostar 
formation. 

Problems such as the p.d.f. and the power spectrum of 
the density field, or the clumps density and mass distribu- 
tion, and the role of the magnetic field will be tackled by 
detailed studies of three-dimensional supersonic MHD nu- 
merical turbulence (Nordlund & Padoan, in preparation). 

It is worth noting that if the present mode of proto- 
star formation predicts very low efficiencies it is likely that 
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APPENDIX: List of symbols 





exponent of the clumps mass distribution 


Cs 


isothermal ambient sound speed 


e 


protostar formation efficiency 


L 


size of the region with r.m.s. velocity a Vi o 


m 


protostar mass relative to Mj tC i 


m c i 


clump mass relative to M r j 


M cl 


mass of the typical postshock clump 


Mj, c i 


Jeans' mass in the typical clump 


Mj 


local Jeans' mass 


M rad 


minimum mass scale for isothermal shocks 


M. 


total mass in protostars 


M tot 


total mass 


M 


Mach number 


M t 


Mach number of the ambient r.m.s. turbulent velocity 


M(100pc) 


Mach number on the scale of lOOpc 
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n 


local number density 


n 


ambient average number density 


n c i 


number density of the typical fragments of mass M r j 


P(m) 


Jeans' mass distribution 


P(x) 


density fluctuations distribution 


P'(m) 


protostar mass distribution 


Pci(mci) 


clumps mass distribution 


tj v 


r.m.s. turbulent velocity 


<?v,0 


ambient r.m.s. turbulent velocity 


P 


volume density 


Po 


ambient mean volume density 


Ttr 


scale of the transition to subsonic turbulent motions 


To 


ambient average temperature 


7~cool 


postshock gas cooling time 


Tdyn(r) 


dynamical time on the scale r 


7~dyn,0 


dynamical time on the scale Lo 


X 


local number density relative to no 



